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Abstract 

Standard lattice-space formulations of quartic self-coupled Euclid- 
ean scalar quantum fields become trivial in the continuum limit for 
sufficiently high space-time dimensions, and in particular the moment 
generating functional for space-time smeared fields becomes a Gaus- 
sian appropriate to that of a (possibly generalized) free field. For 
sharp-time fields this fact implies that the ground-state expectation 
functional also becomes Gaussian in the continuum limit. To overcome 
these consequences of the central limit theorem, an auxiliary, nonclas- 
sical potential is appended to the original lattice form of the model 
and parameters are tuned so that a generalized Poisson field distribu- 
tion emerges in the continuum limit for the ground-state probability 
distribution. As a consequence, the sharp-time expectation functional 
is infinitely divisible, but the Hamiltonian operator is such, in the gen- 
eral case, that the generating functional for the space-time smeared 
field is not infinitely divisible in Minkowski space. This feature per- 
mits the models in question to escape a manifestly trivial scattering 
matrix imposed on all infinitely divisible covariant Minkowski fields. 
Two sequentially related proposals for an alternative lattice formu- 
lation of interacting covariant models in four and more space-time 
dimensions are analyzed in some detail. 
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1 Introduction and Overview 



One of the more useful techniques to generate a Euchdean scalar quantum 
field theory is to take the continuum and infinite volume limit of a model 
theory defined on a finite space-time hypercubic lattice. Let the expression 

schematically represent the generating functional for such a lattice theory. 
Here {hk} is an external field, k = {ki, k2, ■ ■ ■ , kn), kj E Z = {0, ±1, ±2, . . .}, 
labels a lattice site, a denotes the lattice spacing, and n denotes the space- 
time dimension. The expression W(0) denotes the lattice action, which for 
a (p^ model, for example, is conventionally taken to be 

W{(f)) = |F(a)S(0fc. - 0fc)2a"-2 + lml{a)J:<j)la^ + go{a)J:<j)la^ , 

where k* denotes one of half of the nearest neighbors to k and the summa- 
tions run over the sites and nearest neighbors as needed on a finite size lattice. 
The continuum limit includes an increase in the number of lattice sites as 
well as sending a —>■ 0. To facilitate convergence to potentially interesting 
results, cutoff-dependent coefficients {Y, mo, and go) have been included. For 
n = 2, 3 the resultant theory that emerges in the continuum limit is nontrivial 
(non-Gaussian) and is in accord with that computed from a (renormalized) 
perturbation theory in the coupling constant go 0. When n > 5, and pos- 
sibly for n = 4 as well, the theory that emerges in the continuum limit is 
trivial (Gaussian), to wit 

h{x)<P(x) d"x'^ = e^/ h(x)C(x-y)h(y)d"xd"y ^ 

which corresponds to a Minkowski theory that is a (possibly generalized) free 
field [0. The manifest non-Gaussian character of the model on the lattice 
is subsumed into a possible influence on the correlation function C in the 
Gaussian limiting behavior. 

The limiting behavior described above for n > 5, and possibly n = 4 
as well, is characteristic of the general behavior that goes under the name 
of the Central Limit Theorem (CLT)[0]. Although the simplest illustrative 
examples of a CLT behavior are for independent identically distributed ran- 
dom variables, neither independence nor identity of the various distributions 
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are essential for its validity. If we accept for the further sake of argument 
that a CLT mechanism is behind the Gaussian limiting behavior for n > 4 
(including n = 4 for convenience), then we are also drawn to consider Pois- 
son distributions as they are the very kind that resist the vise grip of the 
CLT and lead to non-Gaussian and thereby potentially nontrivial continuum 
limits. By Poisson distributions we mean all those that go under the name 
Poisson, compound Poisson, or generalized Poisson distributions 0, which 
are described below as needed. It is the latter of these distributions that will 
be of physical interest in the present work. 

Once the possibility of Poisson distributions is raised, there are (at least) 
two ways to proceed. From the point of view of Euclidean quantum field 
theory — where we presently find ourselves — a natural procedure, which we 
shall initially discuss, suggests itself, but we shall find it wanting in that al- 
though the Euclidean theory may be non-Gaussian the associated Minkowski 
theory is forced to have trivial scattering 0. Our second construction of a 
Poisson distribution is qualitatively different than the first one, and nontriv- 
ial scattering in the Minkowski theory is not manifestly excluded. Although, 
in hindsight, the second approach may appear more natural, both stand- 
points are included for completeness. We now outline the two alternative 
approaches in more detail. 

One may first be tempted to insist that the entire space-time dependent 
field distribution should be made into a generalized Poisson distribution in 
Euclidean space. This change of distribution is presumed to arise by the ad- 
dition of an auxiliary, nonclasssical [0(/i^)] potential P to the conventional 
lattice action. Assuming for the moment we could find such a P let us ex- 
amine its consequences for potential nontriviality. In particular, we suppose 
that both the Euclidean and Minkowski fields have the property of infinite 
divisibility 010]. For the Minkowski field this means that the generating 
functional for the smeared field operator 



where |Oij) is a normalized state and (pFj{h) is a field operator for each R. As 




has a form implicitly given by 



(0|e*^('^)|0) = hm exp{R{On\ [e^'^^W - 1]\0r)}, 
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a consequence of this structure it follows, for example, that 



(0|y.(/iM/2M/2M/i)|0r>0. 
Combined with Schwarz's inequality we have 

< \{oH9iM92Mf2Mfi)\of\' 

< {(^\^i9iM92M92M9i)\0f 

x{o\^ihMf2Mf2Mfi)\of , 

which for the asymptotic states of the Haag-Ruelle scattering theory implies 
that 

< \{0\(Pouti9l)^outi92)^inif2)^tnifl)\0)'^\'^ 
< {0\Vout{9l)Vout{92)^out{92)Vout{9l)\0)^ 

x{0\^inifl)^in{f2)^inif2)Vinifl)\0f = . 

The vanishing of this expression implies an absence of two-particle elastic 
scattering. Similar arguments may be extended to multi-particle scattering 
including particle production and annihilation 0. Accepting the hypothesis 
that all asymptotic states are covered in this way leads to a unit S'-matrix. 
Thus, insisting that the entire space-time dependent field distribution be a 
generalized Poisson distribution leads to trivial scattering, and so this par- 
ticular hypothesis loses its appeal. However, there is also another route to 
follow. 

Let us assume for the sake of argument that the sharp-time field is well 
defined (a discussion of lifting this assumption appears in the next Section). 
We observe that if the space-time smeared field has a Gaussian generating 
functional then this property is inherited by the sharp-time field as well. If 

we set 

h{x) = fi'(x)5(t) , xGR% s = n- l , 

then, for n > 4, 

At sharp time, e.g., t = 0, the Euclidean theory expectation coincides with 
the Minkowski theory expectation, and thus it follows that 
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where wc have changed g to ig for later convenience. We may diagonahze 
the field at t = so that 

where we have adopted a formal notation for a straightforward Gaussian 
functional integral. The expression ^g{0} represents the ground state func- 
tional, a Gaussian eigenstate of the Hamiltonian operator. 

The Gaussian character of the continuum ground state contrasts with 
the non-Gaussian ground state of the lattice Hamiltonian operator prior to 
taking the continuum limit. The lattice Hamiltonian may be readily deduced 
from the lattice action, and it has the form 

Ha = -|y-^(a)a-^E'— 2 +|y(a)E'(0fc*-0fc)V-2 

+^ml(a)T.'(l)la''' + go{a)T.' (fyla"" + const. , 

where we assume that go(a) > and that the prime on various sums indi- 
cates that only lattice sites with the same "time" value (equal ki, say) are 
included. The "const." is adjusted so that the ground state ^'a(0) of the 
lattice Hamiltonian satisfies 

In view of the nonvanishing quartic interaction, ^a(0) can not be a Gaussian, 
and as a consequence the characteristic functional 

/ e'^''^'=^'='^^|*„(0)|=^H'ci0fe 

is non-Gaussian although, as noted above, the continuum limit of this lattice 

characteristic functional is Gaussian. Thus, wc may reasonably conclude that 
the distribution of field values at sharp time (e.g., t = 0) succumbs to the 
CLT, and in the continuum limit passes from a non-Gaussian to a Gaussian 
distribution. 

Our strategy at this point is to add an (alternative) auxiliary, nonclas- 
sical [0(fi^)] potential P to the lattice Hamiltonian which has the effect of 
changing the ground state in such a way that in the continuum limit the 
characteristic functional for the sharp-time field corresponds to a generalized 
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Poisson distribution. In symbols, in the first approach we adopt, we modify 
the lattice Hamiltonian so that its ground-state characteristic functional has 
the form given by 

J e^^''t»'3,a^\^,{<f))\^U'd(j)k = exp{-J[l - cos(SVfc^7,a^)]|V'a(0)|'n'#fe} , 
where 

I |^„(0)|2n'#, = oo, 

and insist, moreover, in the continuum limit and written formally, that 

= exp(- |{l-cos[0((7)]}|^{0}|2 50), 

where 

J |^{<^}p5</) = oo. 

In writing these expressions we have made explicit use of the expected sym- 
metry ipa{—(p) = i'ai'P) appropriate to the ground state in order to introduce 
the "cosine". The expression in the exponent in both of the last equations 
is often termed the second characteristic it also serves as the genera- 
tor of the truncated equal-time correlation functions, when such moments 
exist. We suppose that the characteristic functional — and thus the sec- 
ond characteristic — is well defined for a wide class of test functions g, say 
g E iS(R'*), the Schwartz space of functions of rapid decrease. How- 
ever, as already indicated, we will insist that / \ipa{(p)\'^'n.'d(l)k = cxd as well 
as = oo, conditions which characterize a generalized Poisson 

distribution for the lattice and in the continuum limit, respectively, as this 
divergence is necessary and sufficient in the present case for the smeared 
sharp-time field operator (p{g), on the lattice and in the continuum, to have 
a purely absolutely continuous spectrum. The contrary situation for which 
/ |'?/'a(0)p n'(i0jt < OO, or especially / iV'l^jP^^ < oo, conditions which de- 
scribe a compound Poisson distribution, will be unsuitable. 

Once having ensured that the sharp-time field has a generalized Poisson 
(and non-Gaussian) distribution, we choose the associated functional \l/p{0} 
to be the ground state of the continuum Hamiltonian and, based on that 
choice, attempt to build multi-time correlation functions in a conventional 
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manner. Insofar as the sharp-time correlation functions are non-Gaussian, 
one may reasonably expect the multi-time correlation functions to be non- 
Gaussian as well. Since the present construction invokes an auxiliary interac- 
tion to generate a Poisson distribution only for the sharp-time fields, there is 
no reason, in general, to expect that the relativistic field will have a Poisson 
distribution for all space and time thereby avoiding manifest triviality of the 
^-matrix as discussed above. 

Faced with the difficulty of analj^ically determining ^a(0) to carry out 
this program we introduce a second model designed so that the continuum 
limit of the time zero field, at least for several test cases, is identical to 
that based on the first model. In the second model wc approximate ^a(0) 
directly thereby enabling us to recover the auxiliary nonclassical potential 
P directly. Given the auxiliary potential we propose a specific form of a 
lattice action that wc believe embodies our proposals for an alternative and 
nontrivial quantization of self-coupled scalar fields such as (y?^ for n > 4 (and 
potentially an alternative quantization for n = 2 or 3 as well). 



It is instructive at this point to introduce a "toy" model to more easily visu- 
alize some aspects of the previous discussion. Let us consider the well-known 
Gauchy distribution in one dimension given by its characteristic function 



Here 6 > denotes an arbitrary scale parameter. The spectrum of x is 
absolutely continuous as signified by the probability density 6/[7r(a;^ + 6^)]. 
Observe the nonintegrable singularity in the second characteristic weight 
function, i.e., / dx/{7rx^) = oo. The two weight functions involved are, 
of course, related to one another. In particular. 



1.1 Illustrative Toy Model 
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While analytically simple, this toy model is special because its moments 
do not exist. A closely related family of toy models is determined by the 
characteristic function 



where p > is another parameter at our disposal P] . The conventional modi- 



fied Bessel function Ki{z) ~ ^7i/2ze~^ for z 3> 1, ensuring that all moments 
of the distribution exist. Moreover, since Ki{z) ~ I/2 for < z <^ 1, it fol- 
lows that as p — > these distributions converge to the Cauchy case. For any 
p > observe that for small x the second characteristic weight function is 
effectively identical, i.e., 6/(7rx^), to that of the Cauchy distribution. Thus 
we are led to the notion of a set of different distributions all of which have 
the same singularity structure of the second characteristic weight function 
at small arguments but which differ from each other by their behavior for 
intermediate and large arguments. In the present case such a set includes, 
e.g., all distributions of the form 



where F G C^, F(0) = 1, F > 0, and + x^)] dx < 00. Although 

we can not analytically specify the L} nonnegative weight function G of 
the characteristic function in the general case, we know that it exists. The 
characterization of a distribution by means of its second characteristic will 
be an important tool in our study. 

One aspect of our approach that does not have an analog in this toy 
model is the relation of the second model to the first, namely the fact that 
these two lattice field theories have different lattice starting points but equal 
continuum limits. 

Before closing the discussion of our toy models it is worthwhile illustrating 
our suggested procedure for constructing multi-time correlation functions. 
For simplicity we confine ourselves to the Cauchy model (even though an 
entirely analogous procedure could be carried through for the toy model 
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with finite moments). In particular, if we identify 

bdx 



J 7r( 



(X2 + 62) 

as the ground-state expectation function, then we are led to represent the 
Hamiltonian operator H. by 



1— + i(a;2 + 62)1/2^1 
dx"^ 2 Qx^ 

2^2 - &2 
■2^X2 ^ 2(x2 + 62) 



For any 6 > it is seen that this Hamiltonian has only one normalizable 
eigenstate, the ground state with zero energy; this paucity of normalizable 
eigenstates is related to the absence of moments for this model. It follows 
that 

fQl^^p'Q^-^nt^^pQ^Q^ = - / e"^^"* , dx , 

^ ' ' ^ Tri Va;2 + Vx^ + 

with analogous expressions holding for higher-order multi-time correlations. 
The procedure illustrated here is analogous to that to be followed in our 
discussion of models in Sections 2 and 3. 

On the other hand, for the sake of completeness, we should also illustrate 
an alternative dynamical extension analogous to the first proposal to have 
a Poisson distribution for fields over all space and time. In this alternative 
dynamical extension we modify the second characteristic so that 

{0\e'P'^e-''^^e'P^\0) = exp[- / J-UP'^'e-'^^e'^'' - l)T^dx] , 

TT J \x\ \x\ 

where in the present case 

^2 (92 

h = _i |- il^;! \x\~^ 

2 dx^ 2 Qx^ 

2 5a;2 a;2 

At t = 0, of course, the two dynamical extensions coincide. The Hamilto- 
nian h has a positive, purely absolutely continuous spectrum, and as a con- 
sequence, the asymptotic behavior of both dynamical expressions exhibits 
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clustering as t ±00. Thus on general grounds both dynamical expressions 
may seem viable. However, as argued earlier, one of these extensions leads 
to manifestly trivial scattering in the field theory case. Consequently, our 
further discussion will be aimed at converting the sharp-time field distribu- 
tion into a generalized Poisson distribution and proceeding in analogy with 
the former of the two dynamical toy models. 

1.2 Summary of Proposed Lattice Model for (pf^ 

Here we summarize the results of the next two sections and present in as 
concrete a form as possible our proposal for a lattice model designed to lead 
to a nontrivial continuum limit and to have the possibility of nontrivial scat- 
tering as a relativisitic theory. In the language with which the conventional 
lattice models are formulated, our proposal is very much the same save for 
two aspects. The first difference is the presence of the auxiliary, nonclassical 
potential that enforces a generalized Poisson distribution on the sharp-time 
field distribution, and the second difference is the wholly new form taken by 
the bare parameters of the theory in terms of the renormalized (or finite, 
nonzero) parameters that eventually determine the theory properly. 

The model is characterized by the Euclidean lattice form of the correlation 
generating function as given by 

S{h) = NoJ expp /ifc0fca" - |r(a)E(0fc. - 0fc)'a(-2) - lm^(a)S02a- 
-goia)J:<f)la'^ - SP(0fc,a)a"] Ud(j)k • 
In this expression ||^ 

where 

A(a) = F(a)-ia-2^ 8-1(1 + 0(3 + 0, 
B(a) = Y{a)-'a-^'4-\l + OC(a) , 

C(a) = M"-2(Ma)2^(i-«)/« . 

In these relations > is a free parameter labeling inequivalent quantiza- 
tions; as we shall show, the cases < < 2 and 2 < ,^ are qualitatively 



m<Pl - B(a) 
[<Pl + C{aW 
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different. For the present, as well as for the most part, we shall confine at- 
tention to < ^ < 1. Our analysis indicates that the bare parameters are 
related to renormalized parameters according to the relation 

ml{a) = {MaY-^m^ , 
go{a) = {Ma^'-'g , 
Y{a) = (Ma)'' , 

where M denotes a mass parameter chosen to ensure that engineering di- 
mensions are preserved. Observe that these three bare parameter renormal- 
izations are multiplicative as opposed to subtractive; this consequence is a 
direct reflection of the profound difference that the auxiliary potential intro- 
duces into the integrand. Moreover, all three bare parameters are small when 
Ma ^ 1, i.e., as the continuum limit is approached. This means that in the 
integration variables indicated the essential support of extends roughly to 
0{{Ma)^^), namely to large values. On the other hand, the contribution of 
the auxiliary potential is signiflcant for much smaller values of the integra- 
tion variables, namely when 0^ ~ C(a) = M'^~^{MaY^^^~^''^^ . In brief, each 
integration variable sees two qualitatively different regions, one for small val- 
ues dominated by the auxiliary potential, and one for large values controlled 
by the terms of the usual model. Of course, "small" and "large" are relative 
terms and may be influenced by a rescaling of the variables of integration; 
for example, if we introduce pk = {MaY(j)k, then it follows that pk has an 
essential support of order one or less. But in this case the influence of the 
auxiliary potential occurs for extremely small values of pk- As the contin- 
uum limit is approached the disparate scales for the two effects become even 
more pronounced. Incidentally, the importance of the very small and very 
large integration domains for each variable would make a naive application of 
either renormalization group arguments or Monte Carlo calculations appear 
rather difficult. 

In a formal continuum limit the density of the auxiliary potential term is 
invariably of the form 

for any value of ^ and any space-time dimension n, where c denotes a suitable 
constant. We will present arguments that show that this term should not 
be interpreted as a renormalization counterterm for the nonlinear interaction 
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term, but rather as a necessary renormalization for the kinetic energy term 
with which it shares an identical scahng behavior. (The analog of the cen- 
trifugal potential arising out of the kinetic energy in nonrelativistic quantum 
mechanics is an appropriate one here.) 

Such renormalization counterterms have been encountered in previous 
models studied by the author, and one may expect that some of the gen- 
eral discussion that held true for previous models may well hold true in the 
present case [0]. In particular, we have in mind the implications of such a 
counterterm for: (i) the presence of a hard-core interaction; (ii) the limiting 
behavior as the nonlinear coupling constant is turned off being not the usual 
free theory but rather a pseudo-free theory which retains the essential effects 
of the hard-core potential; and (iii) the existence of a meaningful perturba- 
tion theory in the nonlinear coupling constant, not about the free theory, but 
instead about the pseudo-free theory. 



2 Lattice Regularized Models: First Formu- 
lation 



In this section we shall present the first of two formulations of a class of 
models each of which we believe has the right behavior so that in the contin- 
uum limit a non-Gaussian Euclidean quantum field theory emerges. Based 
on the discussion in the preceding section it is possible that in Minkowski 
space these models exhibit nontrivial scattering as well. Among such models 
are natural candidates for nontrivial theories of ip'^ in space-time dimensions 
n > 5 and perhaps n = 4 as well. We begin this section with some remarks on 
the difference between locality on the lattice and locality in the continuum. 

Lattice and continuum locality 

Let us illustrate the point of this subsection with a very simple example, 
namely one that deals with Gaussian white noise. Suppose we want to find 
a lattice prescription for the Gaussian functional integral given by 
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The natural lattice formulation of this problem is 

in which the manifest locality — even ultralocality — of the continuum result 
is already apparent in the lattice form. However, it is by no means necessary 
that the lattice formulation needs to be local in order for the continuum 
limit to be local. As an example we study the following alternative lattice 
formulation for the same problem. 
Let us consider 

where S is a positive-definite matrix, is independent of the lattice spacing, 
and has the property that S^Sfc;; = 1 for each fixed /. For general S, the 
lattice exponent is certainly not local on the lattice; however, the continuum 
hmit leads to the same result as before. Specifically, 

a^O " J 

An acceptable example for S is given by 

S,;, = Mn^^i sin2(A;,- - - y% 

where M is a normalizing constant. This expression has the apparent feature 
of being nonlocal. However, in the continuum any two points a finite distance 
apart are actually an infinite number of (continuum) lattice steps apart, and 
so there is no real nonlocality after all. Insofar as the continuum limit is 
concerned, the matrix S^.; is just as good as the choice S^^i. One matrix may 
be more natural than the other, but it would be incorrect to say that one is 
"right" and the other is "wrong" . Since the continuum limit is the ultimate 
goal they are both equally "right" . 

The purpose of this comment is to assuage any concern about continuum 
nonlocality that may be raised in the lattice formulation to follow. What 
appears to be nonlocal is in fact perfectly local in the continuum limit. 
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2.1 Determining the Second Characteristic 

As a preliminary we choose a "time" direction in the Euchdean lattice; this 
will become the time direction when one changes from a Euclidean to a 
Minkowski theory. For this purpose we single out the first component ki of 
the lattice label as the future time direction, but for convenience we shall 
simply refer to ki as the time direction. Evidently there are s ~ n — 1 
space components, k2, ■ ■ ■ ■ kn- remaining. We assume the original lattice is a 
hypercube with an odd number, 2L+ 1, of lattice sites in each direction. Thus 
the total number of lattice points is iV = (2L + 1)", while the total number 
of lattice points in a spatial slice (at constant ki) is given by N' = {2L+ ly. 
We adopt periodic boundary conditions in the spatial directions and Dirichlet 
boundary conditions in the time direction. Let us next introduce a family 
of real constants {f3k} (with k a lattice site label) each of which is positive, 
(3k > 0, and subject to the normalization condition 

KPk = 1 , 

where the sum is over the A^' sites in any spatial slice in the lattice, i.e., only 
n — 1 dimensions of the lattice. In doing the sum we hold ki fixed, and the 
normalization holds for any ki. We shall use a "prime" to signify when a 
sum, or a factor, refers to a lattice of codimension one. A suitable choice of 
constants (3^ is given in the form 

where K > 1 and M' is a normalization factor given by 

M' = [{K - 1)/{K + 1 - 2K-^)Y . 

Armed with these definitions we introduce an ansatz for the weight func- 
tion appropriate to the second characteristic. In particular we let 



where J is a normalization factor (to be chosen below), w{<p), with w(0) = 0, 
is a function designed to control the large field behavior, while the factor in 
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the denominator controls the small field behavior. We choose the parameter 
7 = 1/2 + ^/(2iV'), < ^ < 1, in order, as we shall see, that 

J |Va(0)rn'd0fe = OO, 

while on the other hand, along with a suitable w, 

y'(sWfcari^a(0)rn'd(/.fe<oo 

for all p > 1. In terms of these expressions, the (negative of the) second 
characteristic reads 

Cs{g) ^J' f[l- cos(SV,g,a-)] U',d<j>, . 

Normalization is readily dispensed with. We choose a special {sp) test 
sequence, for example we may choose 

for a suitable constant M. This special sequence is chosen so that 

E'(^f )V = M + 0(a,L-i) 

is finite and nonzero in the continuum limit. We declare that Cs{g^^) = 1, a 
condition that fixes the value of J, in particular, so that 

^ / [1 - cos(EV,g-V)] .2^, Kd^k . 

In order to better ascertain the continuum limit for the second character- 
istic, it is convenient to introduce a new set of integration variables which we 
shall refer to as hyper-extreme spherical coordinates. Based on the definition 

= Kr]k, K > 0, with = K^, let us consider the measure 

= /.^'-^ /.-^'-^ dK' 5{i - T.'4Wi:[(3u-iv^r Kdm 

= 2K-^-^di^dT{'q) , 
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where 

In terms of these variables it follows that 

Cs{g) = 2J2|[1 - cos(«EWa^)] 6'''''^''''^^-^-^ dKdr{ri) , 

and more particularly for the moments that 

X kP-^-^ dK 5(1 - T^'ril) Kdm ■ 

Here we restrict attention to < < 1 which is readily seen to lead to the 
desired conditions that J{T^'4>kgk(i'^y\ipa{4>)\^^'d4>k diverges for p = and is 
finite for p > 1. The case where ^ is larger is treated below. 



First look at acceptable expressions for w{(f)) 

When moments exist, an acceptable expression for is one for which 

moments of the second characteristic are all of the same order of magnitude. 
We can assure this by studying the quotient of any two moments and ensuring 
that it neither vanishes nor diverges in the continuum limit. We focus our 
initial attention on the quotient 

< {T.'gk^kd'Y > ~ /(S' yfc%a-)2/«i-« e-2"'(«'?) dn dr{ri) ' 

Suppose that w — k^/T, then it follows that 

< {J:'gk(l)ka')' > ^ r^ I{^'gkVka')'dr{r]) 
<{T.'gk(t)ka'f> JiJ:'gkVka^ydT{ri) 

{i:',gla^)mgfa^)Jr^lr^fdT{v) 



3cTa' 



ii:',gy)Jr]ldTiv) 



where c denotes a fixed, nonzero, dimensionless constant. Now, as discussed 
below, it is reasonable to assume something like 
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As a consequence it is necessary to choose T oc (Ma)~^* — for some fixed 
parameter M witli tlie dimension of a mass — in order to have a ratio that 
stays finite and away from zero as a — > oo and A^' — > oo. Note that T has the 
mass dimensions of 0^, namely n— 2. It is readily seen that if w = k^/T'^, with 
T oc (Ma)^^*, a suitable behavior is also found. In addition all other moment 
ratios arc maintained correctly with these same choices. Consequently, we 
conclude that the proper choice for w is given by 



for a large class of functions T that may involve dimensional constants, but 
are (effectively) independent of either the lattice spacing a or the number of 
lattice sites in a spatial slice A^'. This remark has strong consequences for the 
form of the regularized coefficients in the lattice Hamiltonian and therefore 
in the lattice action. 

Elementary exercises in many-variable integrations 

There are a number of significant differences in many-dimensional integra- 
tions between Gaussian-like integrands and Poisson-hke integrands of the 
general kind that appear in the second characteristic. For this remark we 
have in mind integrals of the form 



where p > 1. Here, in this and the following exercise section, we do not dis- 
tinguish between A^ and A"', and for simplicity in notation we drop the prime. 
Note that the integrals in question are simple "caricatures" of the more com- 
plicated integrals of interest; nevertheless these integrals will illustrate an 
important point. More relevant integrals are studied below. 

If we use the variables k and {rjk} introduced previously, it follows that 



w{(t)) = T{K-q{May ) 




as well as those of the form 



7g(2p) = 2 j 
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Evaluation of this expression for large N is dominated by the factor k^"^, 
and a steepest descent method can be used to evaluate to a suitable accuracy 
the integral over k. It follows that the stationary point for each integral of 



this form is given to leading order by k = yN/(2A). As a consequence, for 
each different value of A the integrand lives on a different set — a disjoint 
set of field values — as oo \T1\. This fact is responsible for divergences 

in perturbation calculations. For example, if we shift the value of A and 
attempt to calculate 

I'ci2) = 2j ^e^^'^\^~UK5{l - ^vl)^dT]k 

in terms of /g(2p) by means of a perturbation series in AA = A' — A, then 
we are led to the series 

4(2) = Ia{2) - AAIaiA) + |(AA)2/^(6) - |(AA)3/^(8) + • • • 

which exhibits divergences as N ^ oo since Jg(2p)//g(2) oc N^p-'^\ In the 
same spirit if one attempted to evaluate 

by means of a perturbation in B, then one would be led to the series 
I'i = /g(2)-BJg(6) j{^4)da{ii) 

+ \BHam j{^vl?da{r^) + ... , 

where we have introduced 

Here again we encounter divergences, e.g., because Ig{6) / IgC^) while 



and so the net value for the first-order correction relative to the initial term is 
of order A^. In point of fact, the kind of divergences under discussion here are 
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extremely simplified examples — analogous to the "caricatures" spoken about 
above — of the kind of divergences that arise in perturbative calculations in 
relativistic c/?^ theory! 

We now turn our attention to a study of Ip. The expressions that we 
shall study are once again caricatures of the more complicated ones studied 
in previous sections; nevertheless they contain the essential ingredients for 
present purposes. In particular, for p > 1, let us examine 

Ip{2p) = /(E0^re-^^'^^4E0^]-^/^n#, 

= 2 y K^p-ig-^^' dK 5{1 - Er)l) n dr)k . 

Now we see no large power of k in the integrand to lead to A^-dependent 
factors. Indeed, 

/p(4) _ / K^e-^^' dK (5(1 - J:r]l) U dr]k 
Ip{2) ~ Jne-^^' dK6{l - Sr/2)n(i% 

/ K^e"^"^ dK 1 

J Ke"'^'^^ dK A ' 

which is A^-independent and finite. With 

7^(2) Ke-^'"' dK5(l - J:r)l)Udr]k 

a perturbation expansion in AA as before is given by 

7^2) = Ip{2) - AAIpiA) + i(AA)27p(6) - l{AAfIp{8) + ■■■ 
which exhibits no divergences whatsoever. In a similar way 
Ip^^J ^e-^-'-B-'^nt ^^5(1 _ T.4)Tldr)k 
admits an expansion in B of the form 

I'i, = 7p(2) - i37p(6) j\Ent)da{ri) 

+ |i?%(10) j{T.4fda{r,) + .-. . 
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The Ip coefficients are all finite in this case, but we still have 

I i^rjt) da{rj) oc N^' , 
J da{v) ^ . 

Consequently, if we reinterpret B as being proportional to A^, say B — BN, 
then each term of the perturbation series makes perfectly good sense as a 
power series in B. 

These simple examples serve to illustrate one aspect of the profound dif- 
ference that exists between many-variable integrations involving Gaussian- 
like integrands compared with those that arise with Poisson-like integrands. 
However, there is one glaring deficiency in our discussion, namely we have not 
included any examples in which field derivatives play a role. Let us remedy 
that situation. 

Advanced exercises in many-variable integrations 

Let us again return to Gaussian-like integrals, this time some that resemble 
scalar field models, namely 

The subsequent calculation depends on the number of space-time dimensions 
involved. Let us start simply and take n = 1, namely 

Ig^^J e-^«'s(^'=+i-%)'«-'-s«'««;^-i dK5{l - ^vl)^dr)k . 

With B ^l,it follows that k ^ v^, r]l ^ and {rjk+i-rjkf ^ a/N. 
This is the case of Brownian motion and it is interesting to observe that in 
the continuum limit the time derivative of 77 is almost everywhere finite. This 
behavior is unlike that for for which 

< (0fc+i - (pkf >=< n'^iVk+i - Vkf >~ N{a/N) = a . 

It is only when the amplitude factor (k) is included that the paths pass 
from almost everywhere differentiable to nowhere differentiable! Now let 
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us turn our attention to higher dimensions, in particular n > 3 (ignoring 
logarithmic behavior characteristic of n = 2). In that case — and assuming 
for simplicity that A^a" ^ 1— it follows that (f)l ^ a'^^-^), k ^ A^(i-Vn)^ 
Tjf PS! a", and {rj^* — r}kY ~ oT' as well. It follows, for example that 

< — (pkY >~ a~''"'~^\ which shows the divergence of the field difference 
comes about solely from the radius factor (k) and is not at all due to the 
variables (r/fe). 

If we add a quartic interaction term we are led to consider 

= 2 y e-^-'nn,*-V,fa--^-Bn^a--C.^Jlnta^^^N-l ^^^^^ _ ^^2^ ^^7?^ . 

In this case the relative order of magnitude of the terms in the principal part 
of the integrand is given by 

g-A«2s(r?fc.-r)fc)2o"-2-B/t2a"-CK4Sjj*a" ^ ^_AiV-BAr(i-2/")_CAr(2-4/"-) 

which has the interesting corollary that the quartic interaction term becomes 
comparable with the kinetic term when n = 4 and actually dominates the 
kinetic term when n > 5. This is how (strictly) renormalizable and non- 
renormalizable models appear in this abbreviated language. 

Now let us make a comparable study for Poisson-like integrals. Specifi- 
cally, and with p > 1, we consider 

Xp{2p) ^2 J (E/ifcr7fca")2f e-^'^'^(''fc*-^^)'«""'-^'^'""«2p-i d«5(l - Er^^) ndr^^ . 

Our criterion of quality in the present case will relate to the comparable 
magnitude of the even moments. To study this issue we examine 

Jp(4) _ /(S/ifc77fea")^e-^'''^(''**-''^)'''""'-^'^'""/t3rfK5(l - Sr72)nrf% 
Xp(2) ~ /(E/tfe%a")2e-^'''^('?fc*-''fc)'«"-'-^'^''^"«;cifi:(5(l - Hcir/fc ' 

To approximate this quotient, we use the estimates that rjl ^ {rik* — 7]^)'^ 
1/N, which for n > 3 hold just as before. Thus, in a self-evident notation, 
the quotient is approximated by 

Jp(4) _ „ /(Efc,;/i|a"/ifa") < r^gr/f > J e-"^^'""''-^^'"" cIk 
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For small a and large N . the prefactor a"'/N is tiny, and the only way to 
counter this tiny factor is to allow significant weight to large k, values by 
making A and B small. To obtain the desired result for the quotient, first 
imagine that A = 0. In that case B oc 1/N leads to a quotient that is inde- 
pendent of a and A^. If instead B = then we may choose A oc a'^/N which 
again leads to a quotient that is a- and A^-indcpcndcnt. When both terms 
are present the same values prove satisfactory. Observe that no reweighting 
of field support to disjoint sets occurs if either parameter is rescaled by a 
finite amount. This fact is reflected in the statement that an expansion in 
powers of S or of A would not lead to divergences order by order. (If this 
property were to carry over to the type of integrals that arise for the models 
that are the main subject of this paper, then it may well happen that for 
such models a perturbation expansion in an interaction term does not lead to 
divergences order by order.) For completeness, for the very special examples 
of this section, we observe that if we had added another interaction term in 
the exponent, such as C/t'^S^^^a", then the correct choice in that case would 
be C oc a"/A^^. If instead the added term were DK^^Sr/^^a", then the correct 
choice would he D oc a'^^'^/N^'^ ! 

In summary, except for the term involving lattice derivatives, each of 
the factors is correctly given by the observation that k enters in the form 
K^J [Ma)^ /N. For the kinetic term, besides a^""'"^^ / N , an additional factor 
of {MaY is needed as well. Insofar as an acceptable caricature of the lattice 
action suitable for "large" field values in the Poisson case goes, we conclude 
that we may choose the integrand to be 

where and g are cutoff-independent parameters. 
Further analysis of acceptable expressions for 

We note that the examples just studied are not directly relevant for the main 
problems of interest because we have, for clarity, omitted from consideration 
in this simplified model the basic factor Y{'jJ\T,'ij3k-irif]~"' . It is straightforward 
to see that this factor reweights the {rj^} variables already constrained to lie 
on the unit sphere in N' dimensions in such a way that points where t]1 = d^i, 
for each lattice site label I, are significantly enhanced. That is, the measure 
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dr is heavily weighted on the unit sphere along each coordinate axis, and 
significantly reduced away from any coordinate axis. As a consequence, and 
in the absence of field gradients, we have [with ((•)) = J{-)dT{r])/ J dT{r])] 



1 

iV' 



{vWJ - ^^-'y- +i.o.t., 

etc., which crudely represent the fact that only very short range correlations 
exist due to the strong peaking along coordinate axes, and the only factors 
available to characterize that short range dependence are the terms. The 
first equality is an identity, and the second one is scaled to satisfy the first 
as a sum rule. The third equality fails to satisfy a sum rule but is of the 
right order of magnitude. There are other terms not explicitly written and 
denoted by l.o.t. (lower order terms) and these involve potentially long range 
contributions but which are of a very small magnitude, so small in fact that 
they will make no contribution in the continuum limit. In the exercises 
treated above, for which the /5-reweighting factor was omitted, the weight 
of the rj variables was distributed uniformly over the unit sphere, and that 
gave rise to qualitatively different correlation functions than we are here 
discussing. The present correlation functions have the feature that as the 
points separate, i.e., |/c — /| — > L, the correlations effectively vanish, while 
in the uniformly distributed case that was not the case. In the present case 
it is actually essential that the correlations vanish so that our ultimate field 
theory of interest can satisfy the cluster property (here in spatial directions 
only). 

Now let us consider the effect of introducing field derivatives. The effect 
of field derivatives is to introduce additional, and longer-range correlations 
into the direction field. Thus, as a consequence, we may expect that the 
correlations discussed above are transformed into 



1 

e-^"l'^-'l(Ma)^ 
1^ 
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for some mass parameter M . 

We are now in a position to discuss the addition of field derivatives to the 
factor w{(f)) in the second characteristic weight function. Based on the way 
that w enters the lattice potential, it is clear that our original identification 
of = T(0(Ma)*) remains a valid conclusion when nondcrivativc terms 
are involved. The behavior of lattice derivatives for space-time dimensions 
n > 3 is such that the derivative of the direction field r] diverges, even 
without the need for large n values; this conclusion stems from the fact that 
{'nk* — f]kY ~ fJk itself. Thus the classically motivated factor a^""^) is incorrect 
on the lattice (and more particularly in the continuum limit) with respect 
to the factor a~^. For dimensions n > 3, this factor must be renormalized 
away by an additional factor of a^, just as was the case in the "advanced 
exercises". For the nonderivative terms we are instructed to replace by 
0(Ma)*; to augment this recipe in the presence of derivatives we need only 
drop the classically motivated factor in forming the lattice derivative. 
Consequently, even in the presence of derivatives, the renormalization of 
to become 0(Afa)* is the only modification required to obtain a model 
with cutoff-independent parameters! With this prescription accepted then 
the inclusion of terms with derivatives does not qualitatively change the 
discussion given earlier about how to choose w. 

The route to these conclusions has not been straightforward because the 
road traveled is not a familiar one. Indeed, it is quite possible that we have 
been too cavalier in our analysis of the needed renormalization for terms 
involving the field gradient; the appropriate factors can no doubt be deter- 
mined by a more careful study. However, for the sake of further discussion 
we shall adopt our present arguments. 

RelcLxation of sharp time requirement 

Up to this point we have limited the parameter ^ so that < ^ < 1. By 
our consistent use of "cosine" in the second characteristics, it is not difficult 
to see that we could also allow < ^ < 2 with no real change. Beyond this 
value we run into serious difficulty because the lowest order even moment is 
no longer finite. Suppose, however, that we dealt not with the characteristic 
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function for the field but for a "renormalized" cube of the field. By this 
expression we mean to consider 



exp{-j2y[l-cos(EWfc«')] 



}• 



In this relation it is clear that the parameter ^ that enters into 7 in the usual 
way can be extended now so that < ^ < 6 and the indicated expression 
will still be well defined. The larger values of ^ do not support sharp-time 
fields, although they do support sharp-time fields cubed in the manner shown. 
Even larger values of ^ may be entertained by considering still higher-order 
renormalized field powers. In brief, therefore, the class of models can be 
extended beyond those that support sharp-time fields merely by extending 
the range of the parameter C,- 

However, for convenience, we shall continue to restrict our analysis to the 
case < ^ < 1 for which the smeared sharp-time field is an operator. 

Continuum limit 

We now have accumulated enough information in order to take a continuum 
limit for certain of the models in question. First we consider a "primordial" 
model in which w{(f)) = 0, and we consider the second characteristic given in 
this case by the expression 



To make further progress we need to investigate the nature of the support 
of the measure r. Recall that earlier we have asserted that the weighting of 
the measure r is concentrated along each coordinate axis on the unit sphere. 
Therefore, it is not difficult to argue that in the continuum limit, the measure 
dr would be effectively the same as 



where F denotes a suitable normalization. Use of this expression for the 
second characteristic leads to 



C°s{9) = 2 y [1 - cos{KE'7]k9ka')]K 



^ ^ dKdT{r]) . 
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= Fr^T!i{\gi\a'f J[l- cos{v)]v-^-^ dv 





\9i\ 









which in the continuum hmit becomes 



a—*U J 



where {K°)-^ = f \g'P{:>c)\^ d'x. 

Next, as a further example, wc include one particular term in w and again 
take the continuum limit. In particular, we examine the case where 

Csig) = 2J^ J[l - cos{KJ:'r]k9ka')]e-^^''^"'''' K-^-^ dKdT{7]) . 

The same discussion regarding the general behavior of the measure r holds 
again, and we conclude that 

Cs{g) = FJ% J[l- cos{Kgia')]e-^^''^"'^" k'^'^ dK , 

which implies for the continuum limit that 

limCs(^) = K Jd'xJ{l- cos[^^(x)]}e-^^''^' g-'-^ dg , 

where 

= j d'x [{1- cos[^/^'(x)]}e-^''^' g-^-^dg . 

3 Lattice Regularized Models: Second For- 
mulation 

In the preceding section we have made an ansatz for the weight function 
of the second characteristic by choosing the form of il'a{4')- We also made 
considerable progress in characterizing acceptable functions ipa{<P), as prelim- 
inaries to finding the lattice ground state ^a{<P) which are related by formulas 
presented in Section 1. Unfortunately, it is all but impossible to find analytic 
expressions for from the given form of i/ja- Thus to proceed further we 
must find an alternative way to analytically express the lattice ground state 
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SO that we can define the lattice Hamiltonian and the lattice action. We 
will do so by choosing a different starting point, but one which has the same 
continuum limit for the second characteristic as just discussed. 

In our second formulation we do not start with the weight function for 
the second characteristic of the sharp-time characteristic function but rather 
with the weight function for the characteristic function itself. In other words, 
rather than initially assume a form for we shall make an ansatz regarding 

directly. We are unable (and not necessarily interested) to determine 
whether or not the distribution of the sharp-time field is a generalized Pois- 
son distribution on the lattice (as was the case previously), but since the 
continuum limit will agree for the two different starting points it follows that 
the continum limit in the second case is indeed described by a generalized 
Poisson distribution. 

In the present case we postulate the form of the lattice ground state to 

be 

^ nMTcRF^ ' ^^"^ ^ ° ' 

where W includes the necessary constant to ensure that / |^'a(0) pll'^ d4>k = 1- 
This starting point may look fairly standard accept for the presence of the 
product factor in the denominator. Observe in the present case there is 
no appearance of the /? parameters. Our goal now is to show that this 
choice leads to the same continuum limit in the two special cases discussed 
previously. 

First we assume that W{(t)) = c, c a constant, a case that corresponds 
to our "primordial" situation. In that case we consider the characteristic 
function given by 

C {g)-N Je ^^'=[02 + C(a)]a+O/2 ' 

where N°, which arises from c, ensures normalization. By a change of inte- 
gration variables it follows that 

= n'{l - - cos{gg,)] |^^^^4rFT07^> " 
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Observe that 



where F = /(it^ + du is a fixed number depending only on ^ > 0. 

We next choose the a dependence of C(a) so that = 2J^a* for some 
constant J. This leads to 

C(a) =M"-'(Ma)2^(^-^)/^ , 
for some mass M. With this choice for C it follows that 

C'(,) = n'{l - 27-a./[l - cos(«.)] ■ 

Although C(a) may, depending on the value of converge to zero, remain 
constant, or diverge to infinity as a ^ 0, the combination a^*'C(a) oc a^^^^ 
always goes to zero as a — 0. Therefore, we finally conclude that 

limC°(^) =exp(-2J°2 [ {1 - cos[gg{^)]}\g\-'-i dg) 

= exp[-X°y \g{K)\^d'x] . 

If we again assert that C"{g''^) = exp(— 1), then it follows that the continuum 
limit is exactly the same as was previously obtained. Hence the ansatz for the 
second form of regularization has been verified for the choice W — constant. 

Wc next study the case W = {Ma)'^''^Ti'(f)l + c, where c, which is different 
than before, is chosen to ensure normalization. The expression to be studied 
reads 



[0l + C(a)]M/2 



M'^'o'^ dg 



k 

In the present case 

dg 



^^""^ («2 + 1) (1+0/2 ' 
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and in the continuum limit the exponent in the integrand of the last integral 
makes no contribution and thus asymptotically Kc is exactly what it was in 
the case when W — constant! Hence we may determine that the continuum 
limit in the present case is given by 



and once again the value of J is fixed by appealing to the special test sequence 
(7*^. As a consequence, the continuum limit in the present case cooincides 
with the second example studied at the end of Section 2. Observe additionally 
that this equality of the continuum limit has been acheived with exactly the 
same renormalization of the additional term — namely exp[— (Ma)^*E'0^] — as 
was used in the previous section. 

Thus, based on the identical form of the continuum limits and the appear- 
ance of an identical renormalization, we are persuaded in favor of adopting 
this second and alternative regularization based on the choice of the lattice 
ground state ^E'a(0) rather than choosing the weight function ipa{4>) that en- 
tered the second characteristic in the first approach. Although we have not 
studied any further choices of weight functions w and W in the two cases, we 
strongly believe that the second regularization captures the essence of what 
is contained in the first regularization. Hereafter, we choose the second form 
of regularization exclusively. 

3.1 Determination of the Auxiliary Potential 

Having agreed that the ansatz 



provides an acceptable choice for a lattice ground state, we are in position 
to determine the form of the auxiliary potential P. To focus wholly on 
the auxiliary potential let us return to the "primordial" model for which 
= c, namely for a ground state given by 




H',[</.| + C(a)]a+«)/^ 



|\|ol/2 



n'.[0^ + c(a)]a+^)/^- 



29 



With this form of ground state the auxihary potential P may be identified 
by way of the reduced Hamiltonian 



o2 



The result is that 



P(0fc,a) = -^nY[a) a ^^^^^ 

1.2... .-1-2. 4-^(1 + 0(3 + 001 - 2-^(1 + OC(a) 

- ^^"^^"^ [0i + C(a)]2 
^ A(a)0i - B(a) 

- [ct>l-rC{a)Y ' 

which is just the expression previously quoted in Section 1.2. Although we 
have normally assumed units such that = 1, in the preceding formula we 
have incorporated the correct factor explicitly in order to establish — for 
this one time only — its proper role. 

We note in passing that the ground state expectation of P diverges in the 
continuum limit when it formally becomes proportional to 0~^(x). Therefore 
the ground state expectation of the kinetic energy term must also diverge so 
that the sum is zero in the present case. This compensation of divergences 
demonstrates that the auxiliary potential serves as a renormalization coun- 
terterm for the kinetic energy. This same interpretation remains true even 
with the addition of other potential terms. 



3.2 Lattice Model Proposal 

It is now a small step to propose a lattice model for the ip^ theory based on 
the auxiliary potential just identified. Let us first recall the discussion related 
to the toy model in Section 1.1. There we learned that the model without 
moments (based on the Cauchy distribution) led to a one-dimensional po- 
tential analogous to the auxihary potential in the lattice model as identified 
with the help of the "primordial" model {W — constant). We also recall that 
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for the toy model involving the modified Bessel function the weight function 
for the characteristic function has, for very small p values, a functional form 
for small x identical to the Cauchy distribution. This fact suggests that we 
define the model by a lattice action which for small field values is equivalent 
to that of the auxiliary potential and for large field values has the functional 
form of the conventional action up to suitable bare parameters. In particular, 
we choose the generating functional for the lattice Euclidean-space given 

by 

S{h) ^ Noj exp(EMifca" - |r(a)E((/)fe* - cf^kfa^^-''^ - lml{a)J:(l)la^ 

-^o(a)S0ta" - E {[A(a)4 - B(a)]/[4 + C(a)]^}a") n#, , 

where A, B, and C are given above. Observe that we have not chosen to 
define the model by picking W for the simple reason that we do not have any 
analytic proposal for that quantity; rather we have returned to the lattice 
action which preserves its simple analytic form even in the presense of the 
auxiliary potential. Observe also that for small Ma the coefficients of the 
conventional action terms {Y, uIq, and go) are all small and this fact makes for 
a greater isolation of the large and small field behavior. Analogously to the 
case of small p in the second toy model, this isolation of the two regimes gives 
even greater credence to our choice of the auxiliary potential determined from 
the "primordial" model as the form adopted when the conventional terms are 
present. 

It may be natural for the reader to ask why we have not started with the 
second formulation of Section 3 and dispensed with the first formulation of 
the model in Section 2 all together. The reason we have included the first 
formulation is because experience has shown that it is much easier to moti- 
vate unusual forms of the weight function in the second characteristic than 
in the characteristic function itself. Only after gaining familiarity with the 
strange denominators in the weight function of the second characteristic — 
factors which are mandated in order to satisfy the physically dictated moment 
requirements for the weight function of the second characteristic — and seen 
the consequences of these denominator terms in the continuum limit, only 
then do the unusual denominators of the weight function for the character- 
istic function of the sharp-time field itself become reasonable and therefore 
acceptable. 
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A remark on O(A^) symmetric models 



It is interesting to observe that if one deals with an iV-component scalar field 
with an O(A^) symmetry, then the principal change involved is the replace- 
ment of all factors 0| by the 0(iV)-symmetric term SQ,(0fc)^, 1 < a < iV, 
etc., along with a corresponding change of 7 or equivalently ^. How this ad- 
ditional symmetry fits into the present scheme may be partially traced from 
previous work [O . 
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